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The study of the hyperfine (hfs) anomaly in neutron rich nuclei can give a very specific and unique 
way to study the neutron distribution and the clustering structure. We study the sensitivity of the 
hfs anomaly to the clustering effects in the 9 Be isotope using two different nuclear wave functions 
obtained in the three-cluster (a + a + n) model. The results are compared to those obtained for 
9,11 Be in a two-body core + neutron model to examine whether the hfs anomaly is sensitive to a 
halo structure in 11 Be. 
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I. INTRODUCTION 

The study of the hyperfine (hfs) anomaly in neutron 
rich nuclei, in particular, those with loosely bound neu- 
trons can give a very specific and unique way to mea- 
sure the neutron distribution. In a previous paper Q we 
have obtained the values of the hfs anomaly calculated 
in a two-body core + neutron model for the 9,11 Be iso- 
topes. The hfs anomaly e is defined as the sum of the hfs 
anomaly related to the Bohr-Weisskopf effect (esw) 21 
and of the Breit-Rosenthal-Crawford-Shawlow (BRCS) 
correction S Q, e = cbw + ti- 
lt was found in Ref.Q that in Be isotopes, the value 
esw is comparable to the BRCS correction 5. The value 
cbw is very sensitive to the weights of the partial states 
in the ground state wave function and might vary within 
50% depending on those weights. In pj, we found a dif- 
ference of about 25% for the total hfs anomaly value in 
the 9 ' n Be isotopes. 

In the present paper, we calculate the hfs anomaly in 
a three-cluster model of 9 Be and study the sensitivity of 
the hfs anomaly to the clusterization effects. The calcula- 
tions are performed with two different three-body wave 
functions |J, |5| . These wave functions differ from each 
other by the choice of the cluster-cluster interaction po- 
tentials used in their calculations. They are characterized 
by the same partial states but contributing with different 
weights. They both reproduce the 9 Be rms radius and the 
three-cluster (a+a+n) dissociation energy (1.573 MeV). 
We compare the results obtained with these two wave 
functions, and also compare the present results with the 
core + neutron model calculations Q, Q . 

The accuracy of the description of the 9 Be ground state 
wave function can be explored by calculating the mag- 
netic dipole and electric quadrupole moments which are 
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determined by the weights and quantum numbers of the 
states. The 9 Be magnetic dipole moment is independent 
of the radial behavior of the wave function but this is not 
the case for the electric quadrupole moment. Calculation 
of both the magnetic dipole and electric quadrupole mo- 
ments thus provides a precise test of the accuracy of all 
aspects of the ground state wave function (for the exper- 
imental data on the magnetic dipole moment in the Be 
isotopes, we refer to Refs. Q, no experimental data 
exists on the hyperfine anomaly for these nuclei). We 
also compare the results for 9 Be + to that for 11 Be + , to 
investigate the sensitivity of the hfs anomaly to diffuse 
(halo)neutron structures. 



II. THREE-BODY WAVE FUNCTION 



The wave function of the fragment relative mo- 

tion in the three-cluster model (a + a + n) of 9 Be is 
described with the Jacobi coordinates in the method of 
hyperspherical harmonics 0, ll(1 | . It is written as 
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where F^m^ 1 ' 2 is the "active part" of the three-body 
wave function carrying the total orbital angular momen- 
tum L with the projection Mr,. XSMs is the total spin 
function of the whole system with the total spin S and 
projection M$ (here restricted to the neutron spin func- 



tion). The wave function F^/ L xl112 (£i, £2) depends on 
the relative coordinates (hyperradius and relative angles) 
and the neutron spin 
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where l± and I2 are the angular momenta of the relative 
motion corresponding to the £1 and £2 coordinates. The 
sum L = li + I2 gives the total angular momentum; K is 
the hypermomentum. 

The hyperradius p 2 = £f +£f is a collective rotationally 
and permutationally invariant variable, £1 and £2 are the 
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translationally invariant normalized sets of Jacobi coor- 
dinates. We choose the Jacobi coordinates as follows: 
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where r2 3 is the relative coordinate of fragments 2 and 
3, and 1*1(23) is the coordinate of fragment 1 relative to 
the center of mass of the fragments 2 and 3. A23 is the 
reduced mass number for the pair (2,3), similarly, ^.1(23) 
is the reduced mass of the fragment 1 with respect to the 
mass of the subsystem (2, 3). 

The five hyperspherical polar angles are ^5 = {9, x, y} 
where 9 is defined by the relations 
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The choice of three different systems of Jacobi coordi- 
nates leads to three-body wave functions with different 
phase factor, the three different Jacobi coordinates sys- 
tems being connected together by defined rotations. For 
9 Be, the T-basis correspond to choose the neutron as 
fragment 1, the two a-particles as fragments 2 and 3. 
The Y-basis associates the fragment 1 with one of the 
a-particles. 

The values of the hypermomentum are K = l x + l y + 2n 
(n = 1,2,...). The hyperspherical harmonics have the 
form 
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and P"'^ is a Jacobi polynomial. N£ y is the coefficient 
of normalization 
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and the normalization condition for the function ip£ y (9) 
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The charge density distribution of the three-body nu- 
cleus entering the calculation of the electronic wave func- 
tions can be obtained as 
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where po is a normalization factor. In the case of 9 Bc, 
the a-particle density distribution is approximated by a 
sum of Gaussians, with the parameters taken from 
and giving a charge radius equal to 1.676 fm. 

In the proceeding calculations two differing evaluations 
of 0,0 as noted in Section |I] are used. 



III. MAGNETIC HYPERFINE STRUCTURE 

Here, we briefly mention the main points of the formal- 
ism for the hfs anomaly calculations. For more details we 
refer to Q] and references therein. 

The magnetic hyperfine interaction Hamiltonian is de- 
fined by 



H = - J(r) • A(r 



d 3 r, 



(4) 



where J is the nuclear current density and A is the vector 
potential created by the atomic electrons. 

The hyperfine interaction couples the electronic angu- 
lar momentum J and the nuclear one I to a hyperfine 
momentum F = J + I. The magnetic hyperfine splitting 
energy W for a state | IJFMp = F > is defined as the 
matrix element of the Hamiltonian TL. 

The functions F kJ , G kJ entering the matrix element 
are the radial parts of the large and small components 
of the Dirac electronic wave function, with the quantum 
number k = ±(J + |) for J — l e =p h and the orbital 
angular momentum l e . The calculations are preformed 
with a realistic electronic wave function (see 

The magnetic dipole contribution to the hyperfine 
splitting Wuj)ff has the form 
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where ai is defined as 
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with the Z components of the magnetic dipole moment 
M Zi ^ s ^(ri), related to the angular momentum U and spin 
Si of each nucleon; the summation runs over all the nu- 
cleons. 

For an extended nuclear charge we have 
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with D 4 = -yjl [s 1 ® C 2 {h)] and C k q = yj ^Y kq {h). 

The quantity etj can be expressed through the hfs con- 
stant for a point nucleus a^f 1 as 



ai ~ af \l + e B w + 5), 
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Here, fj, =< II | XXffs s « + 9i^%) I H > defines the mag- 
i=i 

netic dipole moment of the point nucleus in nuclear mag- 
neton units /in- The functions Fq J , Gq J are the radial 
parts of the Dirac electronic wave function in the point 
nucleus approximation. 

The hfs anomaly in the Bohr-Weisskopf effect is 



where b 
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The index j is related to the A nucleons; i(l — 3) de- 



and 



notes one of the three-clusters of nucleons. 

(7) 

f/s are the gyromagnetic ratios of the j-th nucleon or- 
bital motion and spin, respectively. 

The hfs anomaly can be approximated by 
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The index i is related to the three-clusters of relative 
coordinate fj and angular momentum li in respect to 
the center of mass system, and with the appropriated 
expressions gf' '><?; and Di, calculated for each cluster. 
The BRCS correction is defined as 5 = 1 - bK a (00) 0. 
This term is defined by the nuclear charge distribution 
and information on the neutron distribution is contained 
solely in esw- 



IV. HFS ANOMALY AND NUCLEAR 
MOMENTS OF 9 BE IN THE CLUSTER MODEL 



To calculate the magnetic dipole moment and the hfs 
anomaly, we use three systems of coordinates related to 
each other by rotation: the T-basis with the neutron as 
cluster i = 1, and the Y(q) systems (q = 1 or 2) with one 
of the a-particles as cluster i = 2 or 3. The index q = 1 
and q = 2 define respectively the rotations 1(23) — > 2(31) 
and 1(23) — > 3(12) between the T and Y-basis. 

The transformation of the hyperspherical harmonic 
function 113 from the T-basis to the Y-basis is defined 



by 
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where (l' x l' y \l x l y ) q are the Raynal-Revai coefficients 



The transformation of the spin part of the wave func- 
tion is written as 



Sl{s 2 S 3 )S x _ ST^ /Ql\Q \1 «'l( s 2 s 3) S i 

X,SM S — / , WxPx/g ASM* 



(13) 
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and 



where 

K 

To calculate the magnetic dipole moment, we consider 
the contribution of each fragment with respect to the 
rest system. In 9 Be, the neutron spin and the a-particle 
orbital motion contribute to the magnetic dipole moment 
as 
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(m s ) , (m y ) q are the expectation values of the spin and 
angular momentum projections of the fragments. (m s ) 
is found in the T-basis and associated with the neutron 
spin; (m y ) q is found in the Y-basis obtained by the ro- 
tation q = 1(2) for i = 2(3) (see JHl and ifTSJll and 
associated with the a-particle orbital momentum. In Eq. 
(|15fl is the weight of the partial state in the T-basis 
for the channel with quantum numbers c = Ll x l y ; and /3|, 
is the weight of the partial state obtained with l|12l) and 
(|13fl for each channel with quantum numbers c' = Ll' x l' y 
in the Y-basis. Here, we mean by channel the set <r of 
quantum numbers characterizing the partial waves con- 
tributing to the ground state wave function. In the case 
of 9 Be represented by the system a + a + n, this set of 
quantum numbers is L, l x , and l y defined in Section ITU 
The nuclear wave function is summed over the hypermo- 
mentum K which is not included explicitly as the matrix 
element does not depend on it. 

= ^gi^ denotes the magnetic 



In Eq. (JT5J), fi 
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dipole moment of the neutron and fiu 

y * 

2 A ^~^ 2 g\ P ^ rn y is the magnetic dipole moment of the a- 
particle in the state with l y as angular momentum in the 
Y-basis. Note that the magnetic dipole moments are ob- 



tained using for the g factors g. 



(n) 



1. 



-3.8260854(90) and 



^4-^2(3) 



The factor | in (|15|l is the center of mass factor 
in the Y-basis corresponding to the a + (a + n) system. 
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We assume that the a-particle angular momentum is 
defined by two orbiting protons, neglecting the spin con- 
tribution of the nucleons. 

The electric quadrupole moment is 



Q = 2]T(//|r l 2 C 2 (f,)|//) ; 



(16) 



where the summation runs over the fragments. 

In the three-cluster model of 9 Be, the hfs anomaly in 
the Bohr-Weisskopf effect is found as the weighted sum 
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where the sign =p corresponds to / = l y ih, and where jj, is 
replaced by the experimental value (fij = — 1.1778(9)//jv 
0) of the 9 Be magnetic dipole moment. 
Here we denote 
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where the ratio takes into account the center of 

mass motion, A, L is the valence fragment mass (the neu- 
tron mass (i = 1) in the T-basis, and the a-particle mass 
(i = 2,3) in the Y-basis). 
So, we can write 
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where is obtained for each channel |^ >. 



RESULTS AND DISCUSSION 



Three-body model of 9 Be 



In the calculations we use two ground state wave func- 
tions of the 9 Be described as the three body a + a + n 
system. These wave functions are obtained with different 
a — a and a — n interaction potentials. 

The first wave function, WF1 Q, is obtained with 
supersymmetric equivalent potential [l4^. With the a- 
particle charge radius 1.676 fm ^JJ this wave function 
gives the value 2.564 fm for the 9 Be charge rms ra- 
dius. This value agrees with the experimental ones, 



TABLE I: The weights of the partial waves 

8 Be-neutron distance r. 
values of the hfs anomaly obtained with WF1 



; , the mean 
the neutron contribution to the 







r n 
fm 


s(n) 
t BW 
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101 


0.498355 


3.274 


-0.0332 


121 


0.315975 


3.159 


-0.0229 


221 


0.152790 


2.976 


0.0123 


123 


0.019661 


3.758 
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223 


0.003674 


3.585 


-0.0286 


143 


0.006608 
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0.001595 
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TABLE II: The same as in Table U calculated with the wave 
function WF2. 
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2.519(12) and 2.50(9) fm (see [TT| and references therein). 
The 9 Be magnetic dipole moment /i9Be = — 1.0531/ijv 
is less by 10% compared to the experimental values, 
-1.177432(3)//jv El and -1.1778(9)^ 0- The cal- 
culated electric quadrupole moment is 53.39 mb which is 
in good agreement with the experimental value 52.88(38) 
mb [H. 

The second wave function, WF2 is obtained with 
the Ali-Bodmer potential |17| and the a-neutron inter- 
action potential (lq. The three-body interaction poten- 
tial is adjusted to fit the three-cluster dissociation energy 
1.573 MeV. The value of the 9 Be charge radius, 2.707, 
fm is overestimated compared to the experimental one. 
The 9 Be magnetic dipole moment pLs Be = — 1.316/xjv is 
10% larger than the experimental value. The electric 
quadrupole moment is Q = 65.42 mb. 

The difference between these wave functions is in the 
radial y dependence of the valence neutron wave function 
(obtained by integration over the x coordinate and the 
angles), the core charge radius and the weights of the 
partial states L, l x , l y for the ground state wave function 
(see Tables |H and HTJ. 

Tables III and ITU show the neutron contribution e^^) to 
the hfs anomaly (first term in (II 7\i ) and the root mean 
square distance of the valence neutron r„ from the 8 Be 
center of mass in each channel <r. 

The contributions to the hfs anomaly in the Bohr- 
Weisskopf effect, and e^M/. and the BRCS correction 
5 calculated with WF1 and WF2 are listed in Table ITTT1 

The contribution of the hfs anomaly from the a- 



5 



TABLE III: The 9 Be charge rms radius, the neutron radial 
distance r n , the contributions to the hfs anomaly from the 
neutron spin (4vV) and from the a-particle orbital motion 
( e< Bw)> t ne values of esw and of 5 calculated for the two wave 
functions WF1,2 (a-particle radius: 1.676 fm) are compared 
to the (core+neutron) results. WF2* refers to another choice 
of the a-particle radius (1.46 fm). 



Value 



core+neutron 



WF1 



WF2 



WF2* 



rms (fm) 
r„ (fm) 

4t (%) 
4t (%) 

e BW (%) 
S 



2.519 



-0.0236 
-0.0451 
-0.0687 



2.564 
3.207 
-0.02281 
0.00085 
-0.02112 
-0.04644 
-0.06756 



2.707 
3.621 
-0.03059 
0.00088 
-0.02882 
-0.04926 
-0.07809 



2.533 



-0.03032 
-0.04634 
-0.07666 



particle orbital motion, tBW (second term in (jl7() ) is 



» 



and its 



small compared to that from the neutron, c B w 
variation is small also (see Table ITTl*)l . 

To explore the sensitivity of the results to the charge 
radius of the a-clusters in 9 Be, we vary the value of the 
a-particle charge radius from 1.676 to 1.636 fm. The last 
value of the charge radius is obtained with regard to the 
negative contribution of the neutron charge distribution, 
(0.34) 2 fm 2 (see [Sill). 

Correspondingly, the 9 Be charge radius changes from 
2.564 to 2.534 fm for WF1 and from 2.707 to 2.678 fm 
for WF2. Owing to the radial behavior of the electronic 
wave functions (see Fig. ^) entering the expressions 10, 
HI U| l. and i|17|) - (|19fl . with a smaller charge radius (rms 
and rmsc) we get a larger value of the hfs anomaly in 
the Bohr-Weisskopf effect. Thus, in our case, we get an 
increase in the hfs anomaly value tsw by 1.6% for WF1 
and 0.5% for WF2. 

The a-a distance is smaller when calculated with WF1, 
so correspondingly, the 9 Be charge radius is smaller and 
the value of the BRCS correction 5 is less (see Table ITTH . 
5 varies within 6% depending on the description of the 
nuclear wave function. 

€bw and the total hfs anomaly e which is the sum 

( e BW + BW + are ma inly determined by the neutron 
distribution. There is no direct correspondence of the hfs 
anomaly value to the nuclear charge radius. For example, 
if we put the radius of the a-particle at 1.46 fm, the 9 Be 
charge radius obtained with WF2 is 2.533 fm and the 
hfs anomaly e=-0.07666% (see Table lilTl WF2*). These 
values are larger than those obtained with WF1 with the 
charge radius 2.564 fm. Thus, even with smaller charge 
radius values one can get larger values of cbw and e. 

Therefore, we can conclude that the hfs anomaly is 
more sensitive to the neutron spatial distribution than 
to the charge distribution of the whole nucleus, and the 
value r n is a crucial parameter for the hfs value in the 
Bohr-Weisskopf effect. 

With the different nuclear wave functions, the r n value 
varies of about 10% and the hfs anomaly £bw of about 




FIG. 1: The electronic wave functions obtained for the dis- 
tributed nuclear charge, G (dashed line), and for the point 
nuclear charge, Go (solid line). 



30%. As the value of the hfs anomaly esw is about two 
times smaller than 5, the total hfs anomaly e varies within 
14% with the choice of the wave function. 

We should also mention that the value of the hfs 
anomaly is very sensitive to the contribution of the dif- 
ferent partial states in the 9 Be ground state wave func- 

(n) 

tion. In particular, the hfs anomaly £ B w m ^ ne channels 
<^ = 121 and q = 221 (that is of 50% of the nuclear wave 
function) is twice as small as 4w m ^ ne channel q = 101. 
Thus, the relative weights of these states are of major im- 
portance. As found in [2lJ the magnetic dipole moment 
in the a + a + n cluster model is also rather sensitive to 
these weights. In the case of 9 Be, three channels mostly 
contribute to the magnetic dipole moment and the hfs 
anomaly, q = 101, q = 121, and q — 221. The magnetic 
dipole moment in this case is [2lT | 



-1.857w 
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1.19Lv 
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(21) 



According to [21j the experimental magnetic dipole 
moment value can be reproduced under the condition 
w 22i < 16 %. Otherwise the calculated value is underes- 
timated. 

Let us analyze the correlation between the weights and 
the values of the 9 Be hfs anomaly, magnetic dipole and 
electric quadrupole moments, noting that the analysis is 
model dependent. 

(n) 

To estimate how the hfs anomaly e B w varies with the 
weights of the different states, we consider the contribu- 
tions of these three channels only (so that ui 



2 

101 



- u> 



121 



w 22i + w res = 1) an( i nn d the weights satisfying the 
experimental value of the 9 Be magnetic dipole moment 
(fXi = — 1.177432(3)/ijv). Under this condition, we get 
the hfs anomaly value plotted in Fig. [21 as a function 



of i 



'101 



for the wave functions WF1 and WF2 (solid and 



dashed lines, respectively). 

Similarly one can express the electric quadrupole mo- 
ment through the weights of these dominant states. For 
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FIG. 2: The variation of the Bohr-Weisskopf hfs anomaly 
value from the neutrons with the weight Wioi for WF1 (solid 
line) and WF2 (dashed line)(a-particle radius: 1.676 fm). 
The dots on the lines are the values in agreement with the 
electric quadrupole moment(see text). The weight wf 2 i is in- 
dicated for each point. 



WF1 we get 



J 221 



Q » 3.157w 2 i - 4.553w 2 21 - 4.527w^ 
+ 2wioi^i2i32.859 + 2wioiw 2 2i32.039 + 

+ 2cJi2lW22l7.314 + Qres, 

where Qres — 11.97 mb. 
For WF2 this relation is 
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Q w 4.204w 2 01 _ 6.635w 2 21 



6.662w 



221 



+ 2o;ioiwi2i41.908 + 2cjioiW221 41.610 

+ 2wi2lW22l9.918-l-Qr.es, 



(23) 



where Q 



fi re 



= 13.11 mb. 

and Q res are the contributions of the ne- 



glected channels respectively to the weights, the magnetic 
dipole and electric quadrupole moments. The square dots 
on the lines in Fig. [21 point the hfs values obtained with 
the weights which also satisfy the experimental value of 
the 9 Be electric quadrupole moment 52.88(38) mb. For 
each wave function we get a few points - two for WF1 
or even four for WF2. On this figure we also report the 
corresponding values of the weight u> 221 {u 221 < 16%). 
One can see the ranges of the hfs anomaly values ob- 
tained with WF1 and WF2 respectively Therefore, an 
experimental estimation of the hfs anomaly in the Bohr- 
Weisskopf effect could give values for the weights of the 
partial states in the ground state wave function. 

The figure shows that the hfs anomaly is a very critical 
quantity to test the nuclear wave function, all other pa- 
rameters being equivalently well described, in particular, 
the electronic part. 



2. Comparison with the core+neutron model 

In Ref. 0, the hfs anomaly for 9 Be was calculated 
within the core+neutron model of 9 Be. The 9 Be ground 
state wave function was given by the superposition of 
states 



| 9 Be(3/2-)) = u; 0+ \l 8 Be(0 + )®n P3/2 } 3/2 -) (24) 
+ uj 2+ \ [ 8 Be(2+)(g>n P3/2 } 3/2 -). 



corresponding to the 8 Be core in the ground (0 + ) and 
excited (2 + ) states with the neutron separation energies 
1.665 and 4.705 MeV, respectively. 

With the weights w 2 + = 0.535 and w 2 + = 0.465 ob- 
tained with the spectroscopic factors from Ref. [2^ the 
magnetic dipole moment is /i = —1.0687 (In and esw, 5 
and e have the values reported in Table ITTT1 The €bw is 
close to the values cbw = —0.0249% Q (obtained with 



the weights w 2 + 



J 2 + 



0.5) and esw 



-0.0243% 



from 23] 

The €bw values obtained for the + state in the 
(core+neutron) model (esw — —0.0440%) and for the 
l x = state of WF1 or WF 2 in the (a + a + n) model 
(e BW = -0.0332% or -0.0402%) are relatively close to 
each other. On the contrary, the cbw values for the dif- 
ferent partial states l x = 2 in the three cluster model 
exceed by a few times (see Tables [I] and |nj the value 
obtained for the 2 + state in the (core + neutron) model 
(e BW = -0.0063%). 

The BRCS correction obtained with the two-body 
wave function is —0.0451%, close to that obtained in 
the three-body calculations. Thus in the core+neutron 
model we get e = —0.0687%, to be compared with the 
values e = -0.06756% and e = -0.07809% obtained with 
the three clusters a + a + n wave functions. 

Thus the clustering effect, revealing itself in the set 
of states contributing to the ground state wave function, 
lead to a variation of the hfs value e of less than 2% for 
WF1 and of about 14% for WF2. Compared to the n Be 
nucleus the difference in the value of the hfs anomaly in 
the Be isotopes is about 25%. This value gives us the 
accuracy of the measurements of the hfs anomaly needed 
to study clustering effects in light nuclei. 

This results corroborates the conclusion in Ref. 0, 
that the value of the hfs anomaly reflects the extended 
neutron distribution in 11 Be and might indicate a neutron 
halo, but the difference for the different isotopes is not 
so pronounced as was found in Ref. 0. 



VI. CONCLUSION 

In the present paper, we have calculated the hfs 
anomaly in the 9 Be + ion with the nucleus described in a 
three-cluster model. The ls 2 2s electronic wave function 
is obtained taking into account the charge distribution of 
the clustered (a + a + n)nucleus and the shielding effect 
of two electrons in the Is 2 configuration. 
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The result of the calculations strongly depends on the 
weights of the partial waves contributing to the ground 
state wave function. Together with the magnetic dipole 
and electric quadrupole moments the value of the hfs 
anomaly can be used to study the clustering effects in 
neutron rich light nuclei. 

The total hfs anomaly is the sum of 6 and 6bw- The 
BRCS correction 5 is only determined by the nuclear 
charge distribution and slightly varies from isotope to 
isotope. The value of the BRCS correction is compara- 
ble or larger than the value of cbw- The hfs anomaly 
in n Bc differs from that in 9 Bc by 25%. The cluster- 
ing effect leads to variations of the hfs value within 15%. 
The calculated magnitude and differential change in the 
value of the hfs anomaly is indicative of the experimental 
precision that must be achieved to study the clustering 
effect and the neutron distribution in neutron rich light 



nuclei. 
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